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Abstract 

The Darboux-Dressing Transformations developed in [T] are here applied to con- 
struct soliton solutions for a class of boomeronic-type equations. The vacuum (i.e. 
vanishing) solution and the generic plane wave solution are both dressed to yield one 
soliton solutions. The formulae are specialised to the particularly interesting case of 
the resonant interaction of three waves, a well-known model which is of boomeronic- 
type. For this equation a novel solution which describes three locked dark pulses 
(simulton) is introduced. 
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1 Introduction 



Many phenomena of wave propagation in one dimensional media in various physical con- 
texts are modelled by systems of Partial Differential Equations (PDEs). Among them, 
integrable models are of special interest and are the focus of our investigation. In fact the 
present article follows our work pQ to which we refer the reader for the general background 
and motivations. Our interest in this subject is particularly motivated by the introduction 
and investigation of boomeron solutions in models of interest in nonlinear optics. These 
solutions are soliton solutions whose behavior differs from that of standard solitons (such 
as for the KdV and NLS equations). In fact they have different asymptotic velocities 
at t = +00 and t = —00. They were first introduced long time ago in [2] (see also [3]) 
in connection with the hierarchy of the matrix Korteweg-de Vries equation, and more 
recently boomerons appeared again in geometry [3] and optics [5], [6], [7J, [E], The 
special properties of boomerons stem from the fact that they are solutions of multicompo- 
nent field equations whose integrability requires a matrix formulation of the corresponding 
Lax pair (see below) together with a richer (non commutative) algebraic structure. The 
characteristic property of this structure is that the flows of the attached hierarchy may 
not commute with each others (see [10] and |llj). The technique of constructing soliton 
solutions we adopt here goes via the Darboux-Dressing Transformations (DDT) which 
yield a new (dressed) solution from a given (naked or seed) solution of the system of wave 
equations of interest. This method may be formulated as a change of the Lax pair of lin- 
ear equations via a transformation which adds one pole to the dependence of its solution 
on the complex spectral variable. Though the formal setting of this construction is well 
known ([12j-[15]). the actual construction is detailed in pQ where the explicit algorithmic 
procedure is given. In particular, there we show that the position of the pole, whether it 
lies off or on the real axis, leads to different expressions of the dressed solution. In its turn, 
as we show here, this difference plays a crucial role in accounting for different boundary 
values of the wave fields, namely boundary values which are appropriate to bright solitons 
(localised pulses with exponential decay) or dark solitons (localised pulses in a non van- 
ishing background). 

To the aim of making this article as much self-contained as possible, in this first section 
we set up the general formalism of matrix evolution equations we have introduced in pQ 
by keeping the same notation as introduced there. In Section 2 the DDT construction, 
which is discussed in details in pQ, is tersely reported. Section 3 is devoted to the soliton 
expressions of both bright and dark type in the general matrix form. However, as for the 
investigation of these expressions, we prefer to focus on one particular reduction. This and 
other reductions lead to models which seem to deserve special attention because of their 
potential applicability. Thus Section 4 is devoted to novel simulton solutions of the well 
know system which describes the interaction of three waves in quadratic media with the 
resonant condition 002 — ^1 = ^3 and k2 — k± = k% on their frequencies and wave numbers, 
respectively. The last Section 5 contains few remarks and conclusions. 
The starting point of our investigation is the following Lax pair of matrix Ordinary Dif- 
ferential Equations (ODEs): 

^ = , iIh = Ti/>, (1) 

where ip, X and T are N x N square matrices, ip = ip(x, t, k) being a common solution of 
the two linear ODEs (pQ), while X = X(x, t, k) and T = T(x, t, k) depend on the coordinate 
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x, the time t and the complex spectral parameter k according to the definitions 

X(x,t,k) = ika + Q(x,t), (2a) 

T{x, t, k) = 2ikC -<rW + a[C, Q[x, t)] , (2b) 

where [A, B] stands for the commutator AB — BA. We notice that, in order to semplify our 
analysis, we have chosen the X and T matrices to be first degree in the spectral variable 
k. This choice has two consequences, namely it lets the two independent variables x and t 
play a similar role so that their interpretation as space and time can change according to 
the physical application, and, secondly, it leads to wave equations which are dispersionless 
when linearized around the vanishing solution. Dispersion terms can be easily introduced 
with only technical changes of the method (see [lj where the dispersive term of Schrodinger 
type has been considered). In the equations (P) a is the diagonal constant matrix 

a = ( 1 N(+)xNM °7V(+)xAT(-) j (3) 
V °JV(-)xJV(+) _1 Af(-)xAf(-) / 

while C is an arbitrary constant block-diagonal matrix, 

C=( C(+) °*%Lf->) , (4) 

where C^ + \ respectively C^~\ are iV"W x N^ + \ respectively x N^~\ constant square 
matrices, with N = + N^~\ and being arbitrary positive integers. The 

potentials Q = Q(x,t) and W = W(x,t) are the off-diagonal and, respectively, diagonal 
block matrices (a superimposed dagger stands for Hermitian conjugation) 



U °AT(-)xJV(-) / ' V iV(-)xJV(+) W< " ) 

where the block U = U(x, t) is a rectangular N^' x matrix while the two blocks 
are square x matrices. Moreover and are, respectively, iV"( + ) x 
and N(~} xJV' - ' diagonal matrices whose diagonal elements , with no loss of generality, 
are just signs, namely 

5( ± )=diag{ S W,..., S W ) }, 4 ±)2 = 1 • (6) 
The compatibility of the two ODEs (pQ) entails the matrix first order differential equations 

U t = U x C(+) - CH U x + U WW + w(-) U (7a) 
= \c {+ \ S {+) u\ , (7b) 

= \c { -\ U S {+) trt , (7c) 

provided the square x matrices and l^M satisfy the conditions 

C (±)t = ^c-Ci)^) (8) 

and 

W (±) t= _ S (±) W (±) S {±) (9) 

The reduced form ([5]) of the matrix Q(x,t) is well motivated by the fact that it captures 
several interesting models of multicomponent wave interactions in weakly nonlinear media, 
see HD1, [TT] and Sections 4 and 5. 
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2 The DDT algorithm 



The method of construction of explicit solutions of the Lax equations ([T]) , and therefore of 
the nonlinear evolution equations ([7]), may be formulated in more than one way |12j-|15j. 
Here we sketch the scheme we developed in pQ where full details are given. 
Let Q<®(x,t) and W^°\x,t) be a pair of matrices with the same block structure (JSJ) of 
Q(x,t) and W(x, t). This assumption defines therefore also the block matrices U^(x,t) 
and (x,t). Let then ip^°\x,t,k) be a corresponding nonsingular (i.e. with nonva- 

nishing determinant) matrix solution of 

v4°) = x(V 0) , ^ (0) = r«V 0) (io) 

where X^(x,t,k) and T^(x,t,k) have the expressions ()2ap and (|2b|) whit Q and W 

replaced by and W<® . Since both compatibility conditions, = V'tcc an< ^ ^ xt = ^Ptx 
are satisfied, (x,t), (x, t) and U(x,t), (x, t) are just two different solutions 

of the same matrix evolution equations ([7|). Consider now the matrix D = ifjt/j^^ or, 
equivalently, the relation 

t/j(x, t, k) = D(x, t, /c)V> (0) (x, t, k) , (11) 

which can be viewed as a transformation if)(°) — > ip. The DDT method consists in searching 
for the dressing matrix D(x, t, k) for a given V> (0) . The expression of the matrix D(x, t, k) is 
found by assigning a priori its explicit dependence on the spectral parameter k. We assume 
that this transformation is near identity and adds one simple pole to the k— dependence 
of the solution t/j^(x,t,k), namely 

D(x,t,k) = l + ^^, (12) 
K — a 

where R(x, t) is the residue matrix and the parameter a is the pole position in the complex 
plane of the spectral variable k. In [1] we derive the expression of the residue matrix R(x, t) 
and we show that this expression crucially depends on the position of the pole k = a, the 
distinction being between the case in which the pole is off the real axis, a/a*, and the 
case in which the pole is on the real axis, a = a* . In the next section we show that this 
distinction is relevant when constructing soliton solutions with vanishing or non vanishing 
asymptotic values of U(x,t) as \x\ — > oo. Here we merely report the following expression 
of the new solution U(x,t), ^W(i,t) in these two cases: 

U = U® - * ,f , (13a) 
£{x, t) 

w (±) = w (0)(±) _ l —— r i . 13b 

£(x,t) 

is the sign matrix Q, is the ATM x AT(±) constant matrix as defined in (jl]) 

with the condition (JSj) and the iV^-dimensional vector z(±) is defined below, see (fT8|) . 
The difference between the two cases is mainly in the expression of the real denominator 
function £(x, t) = £*(x, t). Thus when the pole a is complex with non vanishing imaginary 
part (i.e. a ^ a*), the denominator becomes [Ij 

£(x,t) = £ c (x,t) = * ,J (* (+) ,S (+) z ( +>) - (^.sHj-H)] , (14) 
2{a — a ) 
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while, if the pole a is real (i.e. a = a*), the expression of the denominator is quite different 
from the previous one and reads pQ 



1 f x f 

£ (x, t) = £r(x, t) = £q + - / dx'g(x',t) + / dt'm(xo,t') 

^ Jx J t 



(15) 



where £q , xq and to are arbitrary real constants and the scalar functions g(x, t) and m(x, t) 
are defined by the expressions 

g(x,t) = (zM,SMzM) + (z(-\s(-h^), 

m(x,t) = (z(+\sMc(+hW) - (z(-),5(-)C(-)z(-)). [Lb) 

Moreover, in this second case, i.e. a = a*, the two vectors and are constrained by 
the relation (z^ + \ S^z^) — (z^~\ S^z^) = , which is a consequence of the reduction 
conditions ([5]), j8|) and (|9|) (see [1] for details). The bracket (a, b) in the previous formulae 
stands for the (nonsymmetrical) scalar product of two vectors, 

(a, b) =a^b = Y^ a jbj , (17) 

3 

and the vectors z^ and z^~> are iV^-dimensional and, respectively, TV^-dimensional. 
They are defined by the explicit block expression 



z^\x,t) 



^(x,t,a*) U , (18) 



rS+\x,t) \ _ ^(0)^ + ^ ( 4^ 
.(+) n „A .H 



where the two complex vectors z and z are constant in both cases; they are arbitrary 
if a a* while, if a = a* is real, these constant vectors have to meet the same condition 
satisfied by z^ and z^~\ namely 

(z^, S Mz^)-(zt\s^) = 0. (19) 

As for the expression (|15p . we observe that, because of the conservation law gt = 2m x , the 
differential form gdx + 2mdt is exact and therefore the denominator expression (|15|) can 
be computed as the integral 

I rx,t 

£(x, t) = £ R (x, t)=£ + - / [g(x',t')dx' + 2m(x', t')dt'] (20) 

^ Jxn.tn 



along any curve in the (x, t) plane. 



We end this section noticing that these explicit formulas are meant to serve as the main 
tools to construct soliton- and, by repeated application of DDTs, multisoliton-solutions 
of those systems of multicomponent wave equations which can be obtained by reduction 
of the general matrix PDEs ([7]). However, they have been obtained by algebra and local 
integration of differential equations by leaving out two important ingredients of solutions 
of applicative interest, namely their boundary values and their boundedness. These two 
issues will be addressed in the next sections where explicit expressions of solutions will be 
displayed. 
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3 Soliton solutions 



We construct now explicit solutions of the general matrix equation ([7]) by means of DDT's. 
The general expressions we derive here will be investigated in more details in particular 
cases in the next section. In the present analysis, the matrix-valued function U(x,t), 
which appears in the spectral equation ifi x = Xip , with (j2a|) and (0), is restricted to be 
in one of two different functional classes. The first one is the class B of those (bounded) 
functions U(x,t) which vanish sufficiently fast, for any t, as \x\ — > oo: 

U(x,t) eB if lim U(x,t) = 0. (21) 

Bright solitons are in this class. We point out here that, if U(x,t) is in B, the matri- 
ces W < ^ ± \x,t) which satisfy the equations ([7bj) and ([7c]) are necessarily x-independent 
(possibly t-dependent) when |x| — > oo. In particular we choose 

W {+ \+oo,t) = -fW , W { -\-oo,t) = if U{x,t) e B , (22) 

being two arbitrary constant matrices. We limit our construction to the one soliton 
solution which obtains by choosing in the DDT algorithm the seed solution as the trivial 
solution U(°\x,t) =0,W^ ± \x,t) = =fF^. 

The second class, T>, is the class of those (bounded) functions U(x, t) which do not vanish 
as \x\ — ► oo, as they asymptotically behave like plane waves: 

U(x,t) eV if lim U(x,t) = U± PW (x,t) (23) 

£—►±00 

where the functions U + pw(x, t) and U-pw(x,t) are two generally different plane wave 
solutions of ([7]). If C/ (x, t) is in this class, also the two matrices W^\x, t) which satisfy the 
equations (|7b[) and ([7c]) behave as plane waves, as shown below. The explicit expression 
of plane wave solutions is reported below, but its computation requires some tedious 
algebra which is not reported here. In fact, it turns out that the plane wave solutions 
can be solutions only of two special cases of the general evolution equations ([7]). The first 
(second) case obtains by requiring that the constant matrix (C*~-*) is proportional 
to the identity matrix with the consequence that the matrix W^ + \x,t) (W(~\x, t)) is 
constant, see (|7bl) (see (jTcj) ). However in this case both CW and W^ + \x,t) (CH a nd 
can be easily transformed away so that setting = = reduces the 

system ((7|) to the two coupled equations 

U t = -CH u x + U (24a) 



(24b) 



while, in the second case, namely •* = ' = 0, the system ([7J reduces to the two 
coupled equations 

U t = U x CW + U (25a) 



= U+) , sr(+> c/t sH U 



(25b) 



Moreover, we apply a similarity transformation to the equations (|24p ((|25p) in such a way 
that the nonvanishing matrix (C^) be diagonal (and real, see ([8])) 

C(-) = diagjc^, • • • , c« ,} , 4") = C H* , (26) 
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in the first case, and 



cW=di a g{cS + \...,cW},4 + )=c^, (27) 



in the second case. We also assume, for the sake of simplicity, that c~- ^ if j ^ n and, 
with no loss of generality, that the diagonal part of the matrix (W^) be vanishing, 

namely wt\x, t) = , j = 1, • ■ ■ , iVH (W^\x, t) = , j = 1, ■ ■ ■ , 
In the first case we find that the plane wave matrix solution U{x,t) has the dyadic ex- 
pression 

U(x, t) = exp[i(K^x - n(-h)]A^A(+^ , (28) 
while the x matrix W^~\x,t) is found to be 

W(-\x,t) = -iexp[i(K ( ^x - n ( -H)]B ( -^ exp[-i(K^x - n^-h)} , (29) 

with the following specifications: ' is an arbitrary real constant diagonal N^' x N^~' 
matrix, 

=dmg{k[-\--- , (30) 

with the assumption kj ^ ^ k n ^ if j ^ n; the iV(~ ) x matrix is also a constant, 

real and diagonal matrix 

^-) = diag{^ ) ,.-- ,a$_,} , (31) 
whose entries take the expression 

iVC") / (-) _ (-) \ 

AW and are arbitrary constant vectors of dimension N^ + ' and, respectively, > 
and ' is the i\T(~ ) x off-diagonal matrix whose entries are 

Btf = , B$ = (#» , f 1-1"!) ) > ^ » • (33) 

\ kj kn ) 

Thus these multi-component plane wave solutions are characterized by wave numbers 
kj ^ and frequencies u>j , j = 1, ••■ ,N^~\ and by ./V"( + ) amplitudes and iV^ - ) 

amplitudes Aj , which are the components of the iV"( + )-dimensional vector A^ and, 
respectively, of the N^' -dimensional vector A^~\ We complete the definition of the 
class T> of solutions for the special case (|24p of the general equation (JT]) by adding the 
requirement that the asymptotic plane wave solutions U-pw at x = — oo and U+pw at 
x = +oo, see the definition (|23p . have the same wave numbers kj ^ and frequencies ' 
and may differ from each other only because they have different amplitudes A^ and 
' . This requirement together with the expression (|32p of the frequencies implies that 
the amplitude vectors A^ and a[_ ' at x = — oo should be related to the amplitude 
vectors A+~^ and A^_ ' at x = +oo by the simple relations 

A (+) = j A H = exp(ie(-))^ { T ) , (34) 
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where &( + > is a constant 

AT(+) x AT(+) matrix in the SL(NM) group defined by the 

condition 

$(+)t S (+)$(+) = 5 (35) 

and 0H is a constant real diagonal AfH ) x AfH matrix, 0(~) = diag{$j • • • , 

This family of plane wave solutions (|28p and (j29|) allows for a complete definition of the 

class T> of solutions of the equations (|24p in the first case. 

In the second case the plane wave solutions are described by similar expressions, which 
read 

U(x, t) = A(~U(+)t exp[i(K(+)x - rt+h)} , (36) 

while the iV (+) x AfH 

matrix W {+ \x,t) is found to be 

W( + \x,t) = iexp[-i(K^x - JlWt)]flW exp[i(A"(+)s - fiWt)] , (37) 

with the following specifications: is an arbitrary real constant diagonal A^ + ) x AfH 

matrix, 

=dmg{k{ + \.-- ,k { + ] +) } , (38) 

(+) ^ 

3 

real and diagonal matrix 



v jVt+)- 

with the assumption fcH ^ fcH if j 7^ n; the A^H x A^ + ) matrix is also a constant 



whose entries take the expression 



^ + )=diag{u;i +) ,..., W ( + ( ) +) }, (39) 



nM 1 (+) _ (+) 

^=-4 + Jfc«-(A(-),5(-)AH) £ gW|4^| 2 % C ; +) ) ; (40) 

n=l,n^j V^i ^™ 



AH and AH are again arbitrary constant vectors of dimension AfH and, respectively, 
A/"H and is the A^ + ) x A^ + ) off-diagonal matrix whose entries are 



Bff = , flW = (AH , s(-UH)^)^)4 + > V^TT) ( 41 ) 




c (+) 



!.(+) 
in / 

Again these multi-component plane wave solutions are characterized by A^ + ) wave num- 
bers feH and frequencies , j = 1, • • ■ ,N^ + \ and by AfH amplitudes Am and AfH 
amplitudes A^- , which are the components of the A^ + )-dimensional vector AH and, 
respectively, of the A^ - -* -dimensional vector AH, Also in this case the definition of the 
class V of solutions for the special case (f25j) of the general equation © is completed by 
the requirement that the asymptotic plane wave solutions U-pw at x = — 00 and U + pw 
at x = +00, see the definition (|23p . have the same wave numbers fcH and frequencies 
uij and may differ from each other only because they have different amplitudes A + ) 
and AH . This requirement together with the expression (|40p of the frequencies implies 
that the amplitude vectors A^' and A[_ at x = —00 should be related to the amplitude 
vectors AH a nd A^_ ^ at x = +00 by the relations 

A (-) = $(-)^H j = exp^eH)^^) ; ( 42 ) 
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where ' is a constant matrix in the 5L(iVH) group defined by the 

condition 

$(-)t5H$H = ff(-) , (43) 

and 0W is a constant real diagonal iV"( + ) x iV~( + ) matrix, 0( + ) = diag{9[ + ^ • • • , fljy^}. 
This family of plane wave solutions (|36h and (137p allows for a complete definition of the 
class T> of solutions of the equations (|25p in the second case. 

The main ingredients of the algebraic construction of solutions U(x, t), W^(x, t) are the 
seed solutions (x,t), W WW (x, t) together with the corresponding solution ip^ (x, t, k) 
of the Lax equations (jlOp to which the DDT applies (see (jlip ). In its turn, the DDT 
introduces additional parameters in the new solution, one being the pole a in the complex 
/c-plane, while other parameters come from the residue matrix. In this respect, a word of 
warning is appropriate. A generic choice of these parameters does not necessarily yield a 
solution U(x,t), W^'(x,t) which is acceptable in the sense that i) it is bounded (i.e. it 
is singularity-free in the whole (x, t) plane) and ii) it is localized (i.e. it belongs to the 
same class, either B or T>, of U^°\x, t)). Meeting these acceptability requirements is a 
property of the new solution U(x,t), W^\x,t) that has to be checked, depending on all 
parameters which enter in the construction, including the sign matrices £(+) and SK 
This analysis is quite simple in the scalar case, while in the matrix case, which we are 
presently dealing with, it requires some effort. Here we limit our attention to the case in 
which the seed solution Zj(°\x,t), W^^(x,t) is the simplest solution of the functional 
class to which U^°\x,t) belongs. As a consequence, we obtain explicit expressions of the 
one-soliton solution. Even so, the acceptability requirements mentioned above will be 
investigated only a posteriori and in particularly interesting cases (see section UJ . 

As for the pole a introduced by the DDT, see (fl~2j) . we point out that its position, whether 
on the real axis or off the real axis, depends on the boundary values of U(x,t) which 
characterize the functional classes B and T>. Indeed, in the construction of soliton solutions 
the pole a has a spectral meaning since it belongs to the discrete spectrum of the Lax 
operator L(U) = —iad x + iaQ, see the spectral equation ([1]) with (|2a|) . In general, if 
U(x,t) belongs to the class £>, then the discrete spectrum of the operator L(U) cannot be 
on the real axis and only Darboux matrices D with a ^ a* may give rise to physically 
meaningful solutions. On the other hand, if U(x,t) belongs instead to the class T>, then 
the correspending operator L(U) may have gaps of the continuum spectrum on the real 
axis. In this case, the Darboux transformation, with a = a* lying in one of these gaps, 
yields meaningful solutions. The computation in this last case, which is precisely the one 
of dark solitons, is made difficult by the gap structure of the spectrum which gets more 
and more complicate the higher the two integers AT(+) and AK-) are. 



3.1 Bright solitons 



Let us first consider solutions belonging to the class B. In this case the seed solution of 
([7]) is U(°\x,t) = and W^^^x, t) = ^piK^) = constant matrix. The corresponding 
solution (x, t, k) of the Lax equations (|10p is 



^+\x,t,k) 

v (0)(_) (M,k) 



^(x,t,k)={ * „, f0 ^" ), (44) 
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where the blocks ip^^ (x,t, k) are 

^ {± \x,t,k) = exp[±ikx}exp[(2ikC^ + F^)t] . (45) 

We now follow the DDT algorithm of section [2] by computing the vectors z^(x,t), see 
(fT8l) , whose expression is 

z (±) ( x , t) = exp(±ia*x) (i) 4 ±] , (46) 

where we have introduced the time-dependent matrices 

nW(f) = exp[(2ia*C (±) + F^jt] . (47) 

In the present case, the DDT with a real pole, a = a*, has to be ruled out since it 
leads to a solution U(x,t), W^(x,t) which is not in the class B or V. Thus the pole 
a has to lie off the real axis, a / a*, and the DDT formulas (PJ with E/"( )(a?,t) = 0, 
^ICl^t) = T^ (±) yield the one-soliton solution 

U(x,t) = ~ exp(-2iax)n(-)(t) J B nW t (t)5 (+) , (48) 
W (±) (x,t) = T F (±) - ^ exp(±2fe) [C^\n^\t)B^ ±) U^(t)S^] , (49) 

where 

A(x,t) = (z (+) , 5Wz(+)) - (zH, SH Z H) = 

= exp(2te)(4 +) , n(+)t(t)5(+)n(+)(t)4 +) ) (50) 

- exp (-2te)(4 _) , n(-)t(t)5(-)n(-)(t)zj _) ) . 

Here a = a + ife, namely a and b are real parameters (the real and, respectively, imaginary 
part of the complex pole a). In these expressions we have also introduced for convenience 
the following constant matrices 

5 = 4- ) 4 +)t , B^ = z^. (51) 

The expressions (|48p and are an acceptable, i.e. bounded, solution if, and only if, the 
following inequality 

(4 +) , n(+)t(t)s(+)n(+)(t)4 +) ) , n(-)t(t),s(-)n(-)(t)4- ) ) < o (52) 

holds at any time. For instance, if = 1 and S^' = — 1, this condition is satisfied for 
any vector Zq . With this condition, this multi-parameter solution U(x,t) , W ± {x,t) is 
the one-soliton solution of the system ([7]) which clearly belongs to the class B. Indeed, 
while the expression (j48]) of U(x, t) exponentially decays as \x\ — * oo, in this limit the ma- 
trix W {± \x,t) goes to a matrix whose expression can be easily derived from the formulas 
(1491) together with (|50p . and which depends only on time. This special solution has been 
introduced in [10] for the system ([7]) with an additional dispersive term of Schrodinger- 
type. There it has been shown that these solutions exhibit a rich phenomenology as a 
consequence of the interplay of the matrices C^> and . These soliton solutions may 
well feature processes such as decay, excitation, pair creation and annihilation, and they 
are known, according to their behavior, as boomerons, trappons, simultons (see also section 
4 for additional references). 
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3.2 Dark solitons 



We now turn our attention to the one soliton solution which asymptotically behaves as 
a plane wave and is therefore in the class T>. We limit the present construction of this 
solution to the evolution equations (flHl) , namely to the special case C(+) = = 0. For 

this reason, and to simplify the notation, in this subsection we drop the superscript (— ) 
wherever reasonable, namely we set = C , Cj = Cj , = K , kj = kj , fi(~) = 

Q , ujj ' = and = £?. The construction of the one soliton solution of the equations 
(|25p in the other case = = is quite similar and is not reported. The solution 

obtains by applying the DDT to the plane wave background solution (|28H and (I29p accord- 
ing to the formulae given in section 2, see (|13|) and (|20p . Therefore, we first compute the 
corresponding expression of the solution ip(°\x,t,k) of the Lax equations (|10j) . To this 
aim we more conveniently write the seed matrices Q^(x,t) and W^(x,t) in the form 

Q(°) = GQG -1 , T^ (0) = GWG- 1 , (53) 

where 

G = G(M)=(J eMi{ ° x _ m)] ), (54) 

o s(+)^+)^-)t^(-) \ J 

) j4 (+)i 

This implies that the matrix 



2 =( A("UWt ' W = -U OB' ^ 



$(x,t,k) = G^fattyWfatik) (56) 

satisfies the following pair of ODEs with constant coefficients 

<f> x = iX$, $ t = -iT$. (57) 

Here X and T are N x N constant matrices (N = + iV~( + )) with the following block 
form 

X( ® = [-iA<i-UM -k-K J' (58) 

J W-{ _ iCA (-) A (+)\ _2kc + B-n J ' (5yj 

where all the quantities A^ , B , K , O are defined by the expressions (|28H3ip of the 
naked solution , W^°\ Since the matrices X and T are constant and commute with 
each other, [ X , T] = 0, a fundamental solution of (|57p is 

$(x,t,k) = exp[i(xX -tT)} , (60) 

and therefore the expression of ip(°\x,t, k) reads 

^°\x,t,k) = G(x,t)exp[i(xX -tT)} . (61) 



In order to compute the one soliton solution U , W by applying the DDT ([13]) , the essential 
step is finding the explicit expression of the vectors and z^'. This amounts to 
choosing the constant vectors Zq and Zq , as in (fT8j) , in such a way that the resulting 
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one soliton solution be in the class T>. To this aim we have to find the eigenvalues of the 
matrices X and T as well as their common eigenvectors. Let (ft De an eigenvector of X 
and T and let x an d v be the corresponding eigenvalues, namely 

X<j> = x<t>, T4> = u4>. (62) 

From the block expression (|58p of the matrix X, the eigenvalue equation X <j) = \4> m &y 
be rewritten as 

(k + x + K)</>^ = -i(AW , ^+))A-) 



(63) 



where, <f>^ is the iV^ -dimensional vector defined by our block notation 

</,(■+) 



(64) 



It is plain that both the eigenvector eft = <j>{k), and therefore its blocks (fi( + \k) and cp^~\k), 
and the eigenvalues x = x(k) and v = v{k) depend on the complex spectral variable k. 
Several consequences of these equations easily follow and are reported here below. 

Proposition 1 : if > 1 then there exist N^' — 1 linearly independent eigenvectors 
such that (/>( + ) is orthogonal to the vector A^ + \ i.e (A^ , (f>^) = 0, and cj)(~' = 0, while 
their corresponding eigenvalues, see (|62h . are x = k an d v = with multiplicity — 1. 
Because of this result, the task is that of finding the remaining + 1 eigenvectors and 
eigenvalues. 

Proposition 2 : if </>( + ) is not orthogonal to A^ + \ (A^ , (f>^) ^ 0, then the eigenvector 
takes the expression 

(/>(+) \ / iS^AM \ 

J = V (AM^^A^ik + x + K)- 1 ^ ) ■ (65) 

Of course this expression is not explicit since it contains the unknown eigenvalue x- I n its 
turn, this corresponding eigenvalue x ls one °f the > + 1 solutions of the equation 

iVC-) |i(-),2 

x = k- (#> , sw#)) £ ^uUt^ (66) 

and the corresponding eigenvalue v, see ()62|) . is explicitly given as function of x by the 
expression 

JV(-) ,,(-)|2 

v = _(^(+) , §(+)#)) y g(-) 1 j ' . (67) 

For the sake of simplicity, from now on we assume that the N^' + 1 eigenvalues x which 
solve the equation (|66p and the corresponding N^' + 1 eigenvalues v given by (I67j) . are 
all simple. Moreover we consider below only the case in which the spectral variable k is 
real, k = k*. 

Proposition 3 : if k = k* and the eigenvalue x is real, x = X*> then also the eigenvalue 
v, corresponding to the same eigenvector, is real, v = u*. 
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Proposition 4 : if k = k* and x an d v are eigenvalues corresponding to the same 
eigenvector, then also their complex conjugate x* an d v* are eigenvalues. 



This last result suggests the following notation: N c > is the number of eigenvalues 
X = Xj > j = I? " " " ! ^c> which lie in the upper half-plane, ImXj > 0. We then conclude that 
the N^+N^ eigenvalues of the matrix X(T) are divided as follows: — 1 eigenvalues 
are explicitly given by Proposition 1, 2N C eigenvalues come in complex conjugate pairs 
and the remaining + 1 — 2N C eigenvalues are real. In this respect, and for future 
reference, we observe that, because of the Hermitianity relation 

X* = T,XZ , T 1 " = STS , 

where 



if S = ±1 it follows N c = 0. Moreover, because of the properties (|68p . the following 
orthogonality conditions hold true. 

Proposition 5 : if k = k* , and if <f>j (<f)j),j = 1, • • ■ ,N C , denotes the eigenvector corre- 
sponding to the eigenvalues Xj an d Vj (Xj an d v j)i then 

(fa , Z<t> n ) = , (0 i! S4) = O, j=l,---,N c , n=l,---,N c , (70) 

and 

A,SW = 0, (71) 
The only non vanishing scalar products are therefore (<^- , £</>j) , j = 1, • • • , iV c . 

For future reference, we also observe the following. 

Proposition 6 : if k = k* and \k\ is sufficiently large, then N c = and the + 1 

eigenvalues x(k) = X*(&) are all real and have the following asymptotic value (see (|66p ): 
for one of them x(^) = k + 0(1/ k), while for the others eigenvalues x(k) = —k — 

kj + 0{l/k). 

By decreasing the value of k two real eigenvalues may collide with each other and create 
a pair of complex conjugate eigenvalues. These collisions of eigenvalues give rise to a 
complicate gap structure of the spectrum of the differential operator L(U) = —iad x + iaQ. 
Therefore, for a given real value of k, the generic vector solution of the Lax pair of equations 
(HJ is a superposition of bounded oscillating solutions corresponding to real eigenvalues 
x(k) and u(k), and of exponentially unbounded solutions corresponding to eigenvalues 
x[k) and v{k) which are instead complex. In order to construct a soliton solution in the 
class T>, only a linear combination of these last unbounded solutions can enter the DDT 
formulae of section 2. 

Proposition 7 : if the pole a is real, a* = a, then the vectors and z^ are given by 
the general equation (fT8l) with (|6T|) . namely 



i-\x!t) ) = G(x,t)exp{i(xX-tT)} ( Z J_\ j , (72) 
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where X = X(a) and T = T(a), see (f58|) and ([59j) . The constant vectors Zq and Zq 
are expressed by the linear superposition 

(-) =E(^i+7A-), (73) 

of the eigenvectors of X{a) and T{a) corresponding to complex eigenvalues. Here jj 
and 7j are complex coefficients which are constrained by equation (|19|) . together with 
Proposition 5, to satisfy the condition 

N c 

(4 + \s«4 + >) - (4- ) ,5(-)4- ) ) = Y^hjWj , *h) + c - c -i = • (74) 



As a consequence of this Proposition, the vectors z^' (x, t) and ) (x, t) take the following 
expression 

/ (+)( t\ \ Nc 

{ ) = G MYfri&h + %<&h) > ( ?5 ) 

where we have introduced the complex phases 

Pj =XjX~Vjt , i = V--,N c . (76) 

Moreover, by taking into account the explicit block expression (|65p . we finally obtain the 
expression of the vectors z^ 

N c 

z (+) = i ^( 7i e^ + %eW ) S« A« , (77) 

and 

(78) 

which enter the DDT formulae (|13|) . From these formulae it is also clear that it remains to 
find the denominator function £(x,t) to complete the construction of the soliton solution. 
Since £ x = (1/2)5, see (USD, and g = 2(z ( - + \ S^+'zW) because of the definition (USD and of 
the condition (z^ + \ SW z^) — {zfc~\ S^z^) = , one has to first compute the function 
g(x,t) whose expression reads 

+ E^i 5^i(7i%e*Cft-W + c.c.) + Ef=x(7i7* + «.)] . ^ j 

Note that the RHS of this expression is the sum of exponentials and a term, the last one, 
which does not depend on x and t. The function m(x,t), see (|16[) . is similarly expressed 
as a sum of exponentials and a constant term. However, because of the relations (see 
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(|20|) ) £ x = g/2 , £ t = m, only the constant term in the expression of m(x,t) needs to be 
computed to obtain the explicit expression of £ (x, t) which finally reads 

+ En=x Ef^n=ii^h^- M ~ cc)] + f(A<+) , SMAM) + cc.) 
- t(AM , SMAM) 2 Ef=il7i7* (^ ( ~ } , S(~)C(a + Xj + K)- 2 A^) + cc] . 

(80) 

This expression of £(x, t) greatly simplifies by requiring that it should not vanish for any x 
at any fixed time t. In fact, if we order the eigenvalues Xj = Qj+iPj according to their imag- 
inary part pj as < pi < P2 < • • • < Pn c , the denominator function £(x,t) exponentially 
diverges at both x = ±00 with the leading term — (A^ , A ( - + ^)\'jn c | 2 exp(— 2p^ c x) 

as x — > —00 and 2^-(^ + ^ , A^)\jn c \ 2 exp(2pAr c x) as x — > +00. Since these two 
terms have opposite sign, one has to impose that either 7jv c = or that A f^ c = 0. Iterating 
this argument eventually leads to conclude that a necessary condition that the denomina- 
tor £(x,t) never vanishes (alias that the soliton solution is not singular) is that either all 
coefficients jj or all coefficients 7^ must vanish. Note that this constraint also implies that 
condition (|74p is automatically satisfied. We also observe that these two choices, namely 
jj = and 7j = formally obtain one from the other by changing the sign of all complex 
eigenvalues Xj an d Vj, and therefore we consider only the first choice jj = 0. Thus we 
have the following. 

Proposition 8 : the step-by-step method of dressing a plane wave solution to obtain a 
novel solution amounts to inserting the expressions (see (|77l [7H1 [80]) ) 

N c 

z (+) ^^t^'^+UM , (81a) 
i=i 

n c 

z (-) = , S^A^)e i( - Kx - Q V ^ 7j e^ (a + X j + K)' 1 AH , (81b) 

3=1 

£{x,t) =£ Q - , S(+)AW) Yl Yl lj y* e i{ Pi-P*n) (81c) 

j=l n=l X 3 ~ Xn 

in the general formulae (|13p with 

U^(x,t)=e l{Kx - nt ~>A^A^ , (82a) 
W^ + \x,t) = , W^-\x,t) = - ie mx-Ut) Be -i{K X -m) ( (g2b) 
and the constant matrix = as given by f|33|) . 

We close this section by showing that this novel solution we have constructed is, if locally 
bounded, in the class T>. To this aim we compute the asymptotic behavior of the matrices 
U(x,t) and W^~\x,t) for fixed t and for x —* ±00. This computation is straight and it 
yields the following expressions 

U(X, t) -» e i(Kx-Slt) A {~) A {+)] ; ^ ^ ±QO ^ ( g3a) 

w(-\x, t) -> _ ie i(Kx-nt) B±e -i(Kx-nt) ^ x ^ ±DO ; ( 83b ) 
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where the asymptotic constant vectors A± ^ and matrices B± are given by 

A { ~ ] = AH , A { ~ ] = TA^ , B + = B, B_ = TBT^ . (84) 
Here we have introduced the diagonal unitary matrix T 

r = (a + Xn c + K) (a + X N c +K)- l = diag{e iei , • • • , e i6 »* } , (85) 
where the iV c phases 6j are defined as 

tan(^) = PNc , j = 1, ■ • • , N c . (86) 

a + q Nc + kj 

Therefore, as displayed by our notation, the U and W parts of the solution are asymp- 
totically plane waves of the expected form ex.p[i(Kx — f2i)].A± ^ A^ + ^ and, respectively, 
_^ e i(Kx-nt) g_^_ e -i(Kx-nt) as x -j-QQ w here the asymptotic constant vectors ' and 

A_ ^ and matrices B + and -£>_ are related to each other by the unitary transformation 
(|85p . namely A^_ ' = TA{ \ and, similarly, B- = TB+T^ , where rr^ = 1. In the following 
section we apply the DDT one soliton solution construction to the simplest case = 1, 
iV"(~) = 2, in both classes B and V>. The choice = 1, = 2 leads to the three 

wave resonant interaction model. 



4 The 3WRI model 



The simplest and yet celebrated system of partial differential equations which serves as 
particular example of the results presented in the previous sections is the one describing 
the resonant interaction of three waves in a medium with quadratic nonlinearity. In our 
notation of section 2 this system is obtained by choosing, for instance, iVW = 1 , N^-' = 
2, C« =0, C(-) = C = C* =diag{ Cl ,C2}, Cl + c 2 ,5W = 1,5^ = S =diag{si, s 2 }, 
and by asking that the solution be of the form 

U(x,t) = ( U f^ ) , #)(M) = 0, FT<-)GM) = ( ° S2W*(x,t) \ 

\ U2{X,t) J V SlW[X,t) J 

(87) 

In this case the general matrix PDEs ([7]) reduce to the following system of three coupled 
first order equations (3WRI equations) 

Ulf + C\U\ X = -S 2 W*U2 , u 2 t + c 2 u 2x = S\WU\ , w x = (c 2 - cx)u*u 2 , (88) 

whose integrability was first proved in [16]. It should be pointed out that this system, being 
first order, is covariant with respect to arbitrary linear transformations of the coordinate 
plane (x , t) in itself. This implies that one can give it the form, and attach to the 
independent variables x , t the physical meaning, which is appropriate to the particular 
application of interest. Therefore we discuss the system ([88]) with the understanding that 
all formulae we display here can be easily translated into corresponding results for any 
other form of the 3WRI equations which is covariantly related to it. 

Let us consider first the one soliton solution in the class B. This solution is constructed 
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by mere application of the DDT formulae corresponding to a complex pole a = a + ib, 
namely ([35]) . (H9]l and ([50]) which in this case read 



Abz^ 

U(x,t) = ^— exp (-2iax) U(t)z^ ' 



where 



w(x,t) = f- 4b{C2 A Cl) exp(-26x)[n(t)4- ) 4- )t nt (t)]21 > (go) 

A(x,t) = exp(26x)|4 +) | 2 - exp (-26a;)(4 _) , n^STI^)^) • (91) 

The constant complex parameters which characterize this solution are: Zq , the two- 
dimensional vector Zq ^ and / which comes from the naked solution (see ([49p ) 

Moreover the 2x2 time-dependent matrix TL(t) is (see (|47|) ) 

H(t) = exp[(2ia*C + F)t] (93) 

and the notation [M]2i in the right hand side of (|90p indicates the matrix element M21 
of the matrix M. From the expression (|9ip of the denominator function A(x, i) it follows 
that this soliton solution is singular, due to the zeros of this denominator, if the signs si 
and S2 are s\ = S2 = 1, s\ = —S2 = 1 and s\ = —S2 = —1 and it is instead regular, i.e. 
bounded in the entire (x,t) plane, only if s± = S2 = — 1. In the particular case in which 
the parameter / vanishes, / = 0, this solution has been first introduced in [TB] and fully 
investigated in [TT]-[20|. It describes the interaction of three bright-type pulses, and its 
experimental observation in optics has been recently reported [2T]. The case corresponding 
to / 7^ has been discovered and investigated later in [11] where the soliton behavior has 
been shown to be of boomeronic type [22]. These soliton solutions describe the interaction 
of two bright pulses ui(x,t) and U2(x, t) and a third pulse, w(x,t), which is instead of 
dark type as it does not vanish when x — > ±00. Corresponding to various choices of the 
linear group velocities c% , C2, this solution features a rich phenomenology which has been 
recently shown to be of applicative interest in optical processes in quadratically nonlinear 
media [5]-[8]. The properties of this soliton solution do not need to be reported here since 
they have been detailed in the literature quoted above. 

Let us turn our attention to the soliton solution in the class T>. We consider only those 
solutions which are obtained by applying the DDT technique with a real pole a = a* , 
see (|T21) . and this amounts to dressing a plane wave solution. Thus we apply to the 
3WRI system ([55]) the general formulae displayed in subsection 3.2 as follows. By setting 
= 1 and = A, this being a 2-dimensional constant complex vector, in (I82ap . 
(I82bj) and ([33]), the naked pla ne wave solution takes the expression 

jo) = e mx-^t)] Aj , j = 1, 2 , w<® = -i( ^^-) A 2 ATeW( fe -* 1 >*-( Ua - Wl >t]} , (94) 
■> \k2-k1J 

where k\ and &2 are arbitrary real wave numbers and the corresponding frequencies uj\, 0J2 
are 

uj l = c 1 k 1 + s 2 \A 2 \ 2 (^-^\ , L 02 = C2k 2 + s 1 \Ai\ 2 ( K ^j^ ■ (95) 
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The continuous spectrum of the spectral problem ip^ = X^ip^ , see (jlOp . which cor- 
responds to this plane wave may either coincides with the entire real fc-axis, or may be 
the real axis with one or two gaps (see below). A gap Z g = {k : k_ < k < k + } of the 
continuous spectrum is defined as an interval of the real axis such that, for k £ Z g , the 
3x3 matrix X(k), see (1581) . 





-isiA\ 

X{k) = | -iAx -k-h | (96) 

has one real eigenvalue and two complex (conjugate) eigenvalues, in contrast with the 
continuous spectrum where the three eigenvalues are all real. The very existence of a gap 
depends on the given data k\, k2, A\, A 2 , s\, s 2 . For instance if s± = s 2 = — 1 the matrix 
(|96p is Hermitian and therefore no gap may occur. In the present case the eigenvalues, 
which solve the equation (166j) , are the roots of the third degree polynomial 

V{ X ) = (X~ kf + d 2 ( X - k) 2 + di( X - k) + d (97) 

where 

d 2 (k) = Ak + kx + k 2 , 

dx(k) = 4k 2 + 2k(h + k 2 ) + hk 2 + sM^ 2 + s 2 \A 2 \ 2 , (98) 
do(k) = 2k{si\A l \ 2 + s 2 \A 2 \ 2 ) + sik 2 \Ai\ 2 + s 2 /ci|A 2 | 2 . 

and the eigenvalues may be explicitly given by Cardano formulae. However we skip such 
details and we limit our present discussion to the observation that the spectrum can have 
at most two gaps. This result follows from the expression DIS(/j) of the discriminant of the 
equation V(x) = which comes to be a fourth degree polynomial of the spectral variable 
k, namely 

DlS(fc) = -^=(h - k 2 ) 2 k A + D 3 k 3 + D 2 k 2 + D x k + D , (99) 

where the real coefficients D n depend only on the given data k%, k 2 , A\, A 2 , s\, s 2 . Indeed, 
since the three eigenvalues, namely the roots of V(x): are an rea l if t ne discriminant is 
negative, DlS(fc) < 0, while two of them are complex conjugate if the discriminant is 
positive, DIS(A;) > 0, the three eigenvalues are certainly all real for large enough \k\ 
because the coefficient of the highest power of DlS(fc) is negative (see also Proposition 
6). Therefore, if the four roots of DIS(A;) are all complex (i.e. with a non vanishing 
imaginary part) then there are no gaps in the spectrum. If instead only two roots are real, 
DIS(/c_) = DIS(&+) = 0, A;_ < k + , the spectrum has one gap k_ < k < k + with these 
two roots k± being the endpoints of the gap interval. Finally, if all roots of DIS(A;) are 
real, DlS(k W ) = DIS(^ 1} ) = DIS(A: (2) ) = DIS(^ 2) ) = 0, k W < < k (2) < k ( +\ the 
spectrum has two gaps, k^ < k < k+ and k^ < k < k+ (see Fig. 1). 

In order to proceed further with our construction of the soliton solution, we assume the 
given data k%, k 2 , Ai, A 2 , s±, s 2 are such that the continuous spectrum has at least one 
gap X g and that the real pole a of the DDT lies inside this gap, a£l s . 



We specialize further the expressions (|8Tj) to the present case, namely 

zW = ^, 4 = 7 1, , M . A j , j = 1, 2 , £ = So - HLe- 21 ^ , (100) 
J a + kj + q + tp 2p 
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Figure 1: One or two gaps of the spectrum. In both pictures s\ = s 2 = 1. On the 
left: k\ = —Jv2 = 1, A\ = A2 = 2 and the gap is given by fc_ < k < where 
fc± = ±|a/13 + 16\/2. On the right: fei = 7,k 2 = 2, At = 1, A 2 = 1/2 and the two gaps 
are characterised by fe_ = —4.5, k± = —2.6 and fe_ = —1.5, fej|/ = —0.5. 



with the following specifications: x = Q + ip is the complex eigenvalue of the matrix X(a), 
see ()96p . which lies in the upper half of the complex plane (i.e. p > 0), /3 = — ^ an d f 
is the corresponding T(a) eigenvalue, see (f6"2"j) . given by ([5T|) . namely 

S1C1 l^il 2 ; g2 ? ■ \M> ■ (101) 



a + q + fci + a + q + k 2 + ip 

Inserting these expressions in the general DDT formula (|13|) . and using the notation £q = 
-|7| 2 exp(-2£ )/(2p) and the new variable 

Z = Im{3-t =p(x-Vt)-Z (102) 

finally yields the one soliton expression 



3 ~ ^-J 



Ae <^ jt) !+£__!_ \ i = lj2j (103a) 



This solution is a (iarA; simulton as it describes three pulses travelling on a non vanishing 
asymptotic plateau (\uj\ |A,-| , ki |(c 2 — c\)/{k2 — k\)\Ai\\A2\) and locked toghether 
to travel at the same velocity V, see ()102p . This common velocity is obtained via the 
definition of the variable £, which implies Im v=pV, and it takes the expression (see (llOip ) 

V = SlCl]Al]2 + S2CM2]2 (104) 

{a + q + kxf+p 2 (a + q + k 2 ) 2 + p 2 

This expression is not explicit since it gives V in terms of the eigenvalue \ = q + ip; we 
note however that it implies the following: 

Proposition 9 : if we assume, for definiteness, the ordering c\ < C2, then, according to 
the signs si and s 2 , we find that 

ci < V- < V < V + < c 2 , if si = s 2 = 1 , (105) 
\A I 2 

1/ = ci - (c 2 - 01)7 — ■ ^ 2 — ^ , -00 < V- < V < V+ < ci , if si = 1 , s 2 = -1 , 

(106) 
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V = c 2 + (c 2 - ci) 



\A,\ 2 
(a + q + /ci) 2 +p 2 



c 2 < VL < V < V + < +00 , if si 



V(a) 



1 . 6 




1.58 




1.56 








1.52 






Figure 2: fci = — fe 2 = 1, At = A 2 = 2, ci = 1.5, c 2 = 1.6. On the left: s x = s 2 = 1 and 
ci < F_ < V < V + < c 2 . On the right: si = -s 2 = 1 and -00 < V- < V < V + < c\. 



The results of proposition 9 are a strait consequence of the relation 



1 



+ 



(a + q + ki) 2 + p 2 (a + q + k 2 ) 2 + p 2 ' 
which follows from the imaginary part of the equation 



q + ip = a 



si 



a + q + ki + ip 



\Ai\ 



S-2 



a + q + k 2 +ip 



\A 2 \ 2 



(108) 



(109) 



this being just the present case form of (|66p . We also note that the relation (|108p shows 
again that this soliton solution does not exists if si = s 2 = — 1. 

Figure 2 shows an example in which only one gap is present. If two gaps and T^ 1 



exist, where Tg is characterised by k { ^_' < k < k 



(0 



.(<) 
+ ' 



1 , 2, the results are similar, e.g. 



in the case of s\ = s 2 = 1 the velocity of the simulton V lies in an interval when a varies 
in the gaps; however the function V(a) (together with its range of values) is different in 
the two gaps (see Figure 3). 



V(a) 



-4.5 -4.25 



-3.75 ->2<5 -3.25 -3 -2.75 



V(a) 



1 


5966 




1 


5964 




1 


5962 




-1.2 

1 


59^8 


SiS -0.8 -0.6 




5956 




/ 1 


5954 





Figure 3: k x = 7, k 2 = 2, A x = 1, A 2 = 1/2, c x = 1.5, c 2 = 1.6; si = s 2 = 1. On the 



left: k 



(i) 



-4.5, k [ l ] 



-2.6; on the right: k 



(2) 



-1.5, k 



(2) 



-0.5. In both gaps 



c x < V- <V <V+ < c 2 . 



Finally, it should be pointed out that this solution, as explicitly shown by (|103l) . describes 
the nonlinear interaction of three plane waves whose amplitudes asymptotically, namely 
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at t = ±00, acquire a phase-shift. The expression of these phase-shifts is explicitly given 
by the general formula (|86p which now reads 



tan( 



P 



a + q + kj 



1,2 



(110) 



As an instance, the dependence of the phase-shifts 6\ and #2 on the parameter a inside 
the gap, k- < a < k + , is plotted in Figure 4. 




Figure 4: k\ = — &2 = 1, A\ = A2 = 2. On the left: s\ = S2 = 1. On the right: 
si = — S2 = 1. Here c\ = 1.5, C2 = 1.6 

Similarly, if there exist two gaps Ig and T*jp we observe that the behaviour of the phases 
8\ and 62 as functions of a is different in the two gaps. An example of such dependence 
is displayed in Figure 5. 
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5 Conclusions and remarks 



The Zakharov-Shabat spectral problem, which leads to an important class of integrable 
partial differential equations, f.i. the Nonlinear Schrodinger, Sine-Gordon and modified 
Korteweg-de Vries equations, has been generalized to deal with matrix and vector depen- 
dent variables long time ago. This matrix generalization has been recently revisited in [10] 
and [11] and further generalized in order to derive integrable matrix partial differential 
equations with matrix, rather than scalar, coefficients. This generalization generates new 
integrable systems whose soliton dynamics may have peculiar properties which are known 
as boomeronic and are different from those of standard solitons. In a previous paper pQ 
we have adapted the standard Darboux dressing method to construct soliton solutions for 
a class of boomeronic-type equations. In this paper we have carried out this construction 
and obtained explicit expressions by dressing both the vacuum (i.e. vanishing) solution 
and the generic plane wave solution. We have specialized our formulae to the particularly 
interesting case of the resonant interaction of three waves, a well-known model which is of 
boomeronic-type. For this equation we have constructed a novel solution which describes 
three locked dark pulses. The stability of this solution, with respect to the sign of the 
coupling constants, is not investigated here and remains to be verified. Other potentially 
applicable integrable models are included in the present scheme. One of them, the Dou- 
ble three Wave Resonant Interaction (D3WRI) equation, seems particularly interesting 
in both fluid dynamics and nonlinear optics. In our present scheme, this system of five 
partial differential equations takes the form 

Ult + ClU lx = ~S 2 W*U 2 , U 2t + C 2 U2x = SlWU\ 
VU + C\Vlx = -S 2 W*V 2 , V 2t + C 2 V 2x = SyWV\ 

w x = (c 2 - ci)(it*it2 + sv*v 2 ) 

where s± = ±1 , s 2 = ±1 , s = ±1 are just signs. While these equations are obviously 
integrable, they may also follow, via the multiscale (alias slowly varying amplitude) approx- 
imation, from the resonant interaction of two triads of quasi-plane waves (with a commune 
mode) , whose wave numbers k% , k 2 , q% , q 2 , k 3 and frequencies u\ , uj 2 , v\ , v 2 , u> 3 satisfy 
the resonance condition 

k 2 - h = q 2 - qi = k 3 , u 2 - u x = v 2 - v x = u 3 . 

Regarding the 3WRI and the D3WRI equations, the following remark is appropriate. Since 
they are first order differential equations any nonsingular linear coordinate transformation 
of the (x , t) plane in itself leaves these equations formally covariant. As a consequence, 
in different applications the coordinates x and t may be given different physical meaning. 
Since we have constructed special solutions, this meaning is irrelevant in our present con- 
text. On the contrary, the meaning of these coordinates x and t becomes relevant if one 
addresses the initial value problem, since one has to single out the evolution variable. In- 
deed in our present dressing construction, the two linear equations of the Lax pair (pQ) have 
been treated on the same foot. Instead, the solution of the initial value problem requires 
that one of the two Lax equations plays the role of the spectral problem which provides 
the spectral Fourier-type data, while the other linear equation entails the evolution of 
such data. In this respect we point out that the initial value problem associated to the 
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variable t as evolution variable (time in the usual physical meaning) is still unsolved. This 
is due to the fact that the matrix T(x, t, k), see (|2bp . depends asymptotically, as x — > ±00, 
on the unknown matrix U through the asymptotic values of the matrices see the 

system (|7|). Finding the solution of this open problem is left to future investigation. 
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